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II ■ Abstract. We study real photons produced in heavy ion collisions at RHIC, and we 

Mh| calculate their spectrum and its azimuthal momentum anisotropy. The photons from 

a variety of sources are included, and the interplay and the time-evolution of those 
sources are modelled in a full 3D hydrodynamic simulation. We quantify the V2 of 
thermal photons produced in ideal and viscous fluids, and the consequences of using 
^ ■ different initial conditions are explored. 
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1^ [ 1. Introduction 

o. 

Real and virtual photons are penetrating probes, and thus their measurement can reveal 
much of the internal dynamics of relativistic heavy ion collisions. A complication, 
however, is that photons are emitted throughout the space-time history of the nuclear 
collisions. In addition, several sources of electromagnetic radiation exist, and these may 
be either signal or background, depending on the physics being pursued. Therefore, 
a complete picture needs to address all these known sources, and to consider their 
evolution in a realistic scenario of the nuclear collision. We concentrate here on the real 
photons measured in heavy ion collisions at RHIC: these will include pQCD prompt 
photons, thermal photons from the quark-gluon plasma (QGP) and from the hadron 
gas (HG), and photons from jets interacting with the hot and dense medium. We 
calculate photon momentum spectra and their azimuthal anisotropy, and consider the 
emission from media in and out of equilibrium. 

2. Hydrodynamical evolution and results 

The nucleus- nucleus collision is modelled in this work by using MUSIC [1]: a three- 
dimensional, hydrodynamic simulation of relativistic nuclear collisions. Importantly, 
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MUSIC can be currently modified to include shear viscosity, 77, and fluctuating initial 
conditions (FIC). 

Quite generally, the hydrodynamic evolution traces a space-time locus which 
satisfies the conservation laws for the energy-momentum tensor and the baryon current: 
d^Tf"" = 0, and a^ J^ = 0. In the ideal case, T^^^ = (e + P) ut^u" - Pg'"' and 
Jg jjjg^i = TT-B^i^- Note that P is pressure, e is the local energy density, and n^ is 
the local baryon density; these quantities are related through the equation of state. 
Also, M^ = (7, 7v) is the local flow velocity relative to some fixed reference frame. In 
cases where the dynamics are dissipative, the conservation laws are modified. More 
specifically, T^'^^i]) = T^^eai + 'K^^irj). These details will not be discussed here, but 
can be found in Refs. [21 El H]. However, it is important to state that the viscous 
corrections will affect not only the bulk evolution but also the microscopic distribution 
functions, /q. A popular ansatz [5] for this modification is /o — ^ /o + ^/, where 
^/ = /o (1 =t /o)p'^P'^7r^i//2 (e + P) T^ (± for bosons/fermions) . For processes where 
production yields are amenable to a rate equation formulation, the importance of viscous 
corrections are thus established by using the out-of-equilibrium distribution functions 
above. It is important to realize that the viscous calculations demand a complete 
revaluation of all the underlying photon rates: the details will appear elsewhere [6]. 
The net yield of photons produced in relativistic heavy ion collisions is then obtained by 



lolal ^ 

Ihemial ^ 

jel+QGP 

direct 




Tl/S ^ i/47t - 



Pt (GeV) 



Pt (GeV) 



Figure 1. Left panel: The yield of photons produced in Au-Au collisions at RHIC 
energy, in the 0-20% centrality class. Shown separately are the pQCD direct photons 
from the primordial nucleon-nucleon collisions, calculated at NLO JJ (all scales set to 
Pt/^), the photons from the jet-plasma interaction (SI 1^, and the "thermal photons" 
discussed in the text. Right panel: The net thermal photon spectrum, with and without 
the viscous corrections (here rj/ s = l/47r, and s is the entropy density). 



integrating the rates for photon production with a time-evolution approach. Finally and 
importantly, this model has to be constrained by data, hadronic and electromagnetic. 
In this context, MUSIC has been shown to reproduce the measured characteristics of 
hadronic spectra. The photon rates for equilibrium hadronic ensembles exist for the 
QGP |10] and the HG [11] phases. The real photons generated by an ideal hydrodynamic 
evolution are considered first, and the results are shown on the left panel of Figure [D To 
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set the scale, only a few photon channels are shown and the fragmentation component of 
primordial nucleon-nucleon collisions has been assumed suppressed [12] . The right panel 
quantifies the influence of viscous effects on the net thermal photon spectrum. Even for 
intermediate values of the photon transverse momentum, the changes are non-negligible. 
Going beyond simple one-body spectra, it is useful to state that the azimuthal anisotropy 
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Figure 2. Left panel: The V2 coefEcient for thermal photons is computed as a function 
of transverse momentum, after an ideal hydrodynamical evolution. The separate 
contributions are the V2 from the photons emitted in the QGP phase, in the HG 
phase, and the net V2- Right panet The effect of viscous corrections on the net 
thermal photon V2 is shown. 



of the photon momentum distribution has been shown to be a probe sensitive to the 
nature of the underlying degrees of freedom [131 113] • This fact is made clear in Figure 
El where the left panel shows the value of the f 2 coefficient for thermal photons only. 
Recall that 

V2= r d<p cos (20) (d^N/d^prdy) / T d<p d^N/d^prdy (1) 

The net V2 as a function of px is shown, along with the individual QGP and HG 
coefficients, using averaged initial conditions (AIC, as discussed later). The right panel 
shows the effect of a finite shear viscosity coefficient on these quantities, after the 3-1-1 
dimensional evolution. The influence of viscosity on this observable is remarkable, but it 
unfortunately acts to make a small signal even smaller, at least for conditions prevalent 
at RHIC. Thus, as is the case for hadrons, the momentum asymmetry of real photons 
is reduced by the shear viscous effects. 

The results of the calculations reported on up to now were all done assuming 
smooth (or averaged) initial state conditions in coordinate space. However, a shift 
in paradigm occurred when recent analyses linked non-smooth, or fluctuating, initial 
conditions with the appearance of non-vanishing odd components of the momentum 
azimuthal asymmetry [15]. The photon spectrum has been studied with fluctuating 
initial conditions (FIG) [16] and will not be shown here for the sake of brevity. We 
rather show results for the V2 of real photons, starting from FIG, in Figure [3l In this 
low centrality class, one observes that the FIGs raise the photon V2, as they do for 
hadrons [2]. 
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Figure 3. Left panel: The different V2 coefficients obtained with fluctuating initial 
conditions (FIC). Right panel: The effect of AICs and FICs on the net V2- 



In conclusion, the results presented here constitute the first photon yield and V2 
obtained in a realistic 3D viscous hydrodynamics approach. Our follow-up work will 
include more calculation details and applications to the LHC. Clearly, the continuing 
analyses at RHIC, the advent of LHC photon data, and the new quantitative link to 
realistic viscous hydrodynamics investigated here will fuel theoretical investigations for 
some time to come. 
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